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Abstract 

The exactly solvable BCS Hamiltonian of superconductivity is considered from sev- 
eral viewpoints: Richardson's ansatz, conformal field theory, integrable inhomogenous 
vertex models and Chern-Simons theory. 



1 Introduction 



The exact solution of a simplified version of the BCS model, obtained by Richardson in 
1963, has received in the last few years growing attention due to its physical applications 
to Condensed Matter and Nuclear Physics, and its connection with integrable models, con- 
formal field theory and Chern-Simons theory. In this article we briefly review the subject 
from a historical perspective, focusing on the relationships between various approaches. 

2 The BCS model of Superconductivity (1957) 



In 1957 Bardeen, Cooper and Schrieffer proposed a model to describe the superconducting 
properties of some metals at low temperatures In its simplest form the BCS Hamilto- 
nian is given by 



Hbcs= J2 - gdY,c]+c]_Cf_Cj'+ (1) 

where Cj^± ( resp. cj_|_) is an electron destruction ( resp. creation) operator in the time- 
reversed states I J, lb) with energies ej, d is the mean level spacing and g is the BCS dimen- 
sionless coupling constant. The sums in (||) run over a set of Q doubly degenerate energy 
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levels ej{j = !,...,(]). One assumes that the energy levels are all distinct, i.e. Si ^ ej 
for i j. The Hamiltonian (||) is a simplified version of the reduced BCS Hamiltonian 
where all couplings have been set equal to a single one, namely g, and it describes the pair 
creation and annihilation between electrons belonging to different energy levels. 

In their historical paper BCS gave an ansatz in the Grand Canonical emsemble (g.c.) 
for the ground state of this Hamiltonian which reads 



\BCS) oc exp 



E 



|0) 



(2) 



where Uj and Vj are the variational parameters given by the formulae 




where is the chemical potential, Ej = [{sj — fi)^ + A"^]^^"^ is the quasiparticle energy and 
A is the superconducting gap. The g.c. state @ is asymptotically exact in the limit when 
the number of electrons goes to infinity. 



3 The Projected BCS ansatz (50's-60's) 



Soon after the BCS work, the pairing model was applied in Nuclear Physics, but it became 
clear that a canonical version of the BCS state should be more appropiate to describe 
nuclei with small number of nucleons [^. This gave rise to the so called Proyected BCS 
ansatz (PBCS) which is given by 



\PBCS)^\^^^c[^cl_j |0) (4) 

which has exactly fermion pairs. Again, Uj and vj are variational parameters but the 
equations fixing them are much more complicated than those of the g.c. case @. Of course, 
in the limit when N ^ oo the PBCS results agree with those obtained with the BCS state. 
The PBCS state was applied in the 90's to study ultrasmall metallic grains. A peculiarity 
of the PBCS ansatz is that, having a fixed number of fermions, the superconducting (SC) 
order parameter vanishes identically, i.e. (cj _(_c^ _) = 0, which thus obscures the nature 
of the superconducting correlations. There are however alternative definitions of SC order 
parameters in the canonical emsemble, which converge to the g.c. one for large systems. 
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4 Anderson's conjecture (1959) 



In 1959 Anderson made the conjecture that superconductivity must dissapear for metalhc 
grains where the mean level spacing d, which is inversely proportional to the volume, is of 
the order of the SC gap in the bulk A Q . A simple argument supporting this conjecture 
is that the ratio A/d measures the number of electronic levels involved in the formation 
of Cooper pairs, and thus when A/d < 1 there are not energy levels correlated by the 
pairing interaction. Apart from some theoretical studies this conjecture remained largely 
unexplored until the fabrication of ultrasmall metallic grains. 

5 Ultrasmall metallic grains (1996-97) 

In the years 96-97 Ralph, Black and Tinkham (RBT), in a series of experiments, studied 
the superconducting properties of ultrasmall Aluminium grains at the nanoscale These 
grains have radius ~ 4-5 nm, mean level spacing d ~ 0.45 mev, Debye energy ~ 34 
mev and charging energy Ec ~ 46 mev. Since the bulk gap of Al is A ~ 0.38 mev one 
meets the Anderson's condition for the possible non existence of superconductivity, namely 
d > A. Moreover the large charging energy Ec implies that these grains have a fix number 
of electrons, while the Debye frequency gives an estimation of the number of energy levels 
involved in the pairing, namely 0, = 2ujD/d ~ 150, which is rather small. Among another 
things RBT found a interesting parity effect, similar to the one happening in Nuclear 
Physics, where grains with an even number of electrons display properties associated with 
a SC gap , while the odd grains showed gapless behaviour. 

These experimental findings produced a burst of theoretical activity focused on the 
study of the reduced BCS Hamiltonian (||) with equally spaced levels, i.e. ej = jd. Some of 
the approaches used to study this model have been: i) g.c. BCS ansatz proyected in parity, 
ii) the PBCS ansatz, iii) Lanzcos method up to O = 23 energy levels, iv) Perturbative RG 
methods, v) Density Matrix Renormalization Group (DMRG) up to = 400 levels, etc. ( 
for a review on this topic consult Q). 

6 Richardson's exact solution (1963) 

Surprisingly enough the aforementioned theoretical works were done in the ignorance of the 
existence of an exact solution of the reduced BCS model, obtained in 1963 by Richardson 
and his collaborator Sherman (1964) and further developed and generalized to other 
models in a series of papers in the 60's and 70's by Richardson himself. In order to 
describe Richardson's solution it is convenient to define the hard-core boson operators 
hj = Cj-Cj^j^, bj = _, Nj = bjbj, which satisfy the commutation relations, [bj, = 

The Hamiltonian (||) can then be written as 
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Hbcs = E - 5d E , (5) 



Richardson showed that the eigenstates of this Hamiltonian with N pairs have the (unnor- 
malized) product form 



\N)n = UtiBuW^c), B, = Ef=i2^bl (6) 

where the parameters {i^ = 1, ■ ■ ■ , N) are, in general, complex solutions of the coupled 
algebraic equations 

'ad^ ^ p - p ^ ^ 2r • - p ' ^'''^ 

y" A*=i(^i.) 3=1 

The total energy of these states is given by <?(A^) = ^Yl!v=\ The (normalized) states (^) 
can also be written as IHI 



|iV)^ = ^ E ^%-i,...,J^)&,V--fcL|vac) (8) 



where the sum excludes double occupancy of pair states and the wave function il) takes the 
form 



N , 

^%u---.jn) = Y.\{w: — (9) 

The sum in @ runs over all the permutations, V, of 1,---,A'^. The constant C in (^) 
guarantees the normalization of the state Q. The BCS Hamiltonian can be given a spin 
representation in terms of the operators = 1/2 — Nj,tj' = bj and tj = bj which provide 
a spin 1/2 representation of the SU{2) group associated to each level j. The Hamiltonian 
@ can then be written as 

Hbcs = - E " y(^^ + ^" + ^tes (10) 
i 

where the matrices = X]5^=i {o, = 0, +, — ) satisfy the standard SU{2) algebra, with 
Casimir T • T = T^T^ + \{T+T- + T-T+). 
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7 Integrability of BCS (1997) 



The integrability of the reduced BCS model was established by Cambiaggio, Rivas and 
Saraceno in 1997 M. These authors found a set of operators 



R, = -t^-gd j: (i = l,...,fi) (11) 



which commute among themselves and with the BCS Hamiltonian dlOj), which in fact can 
be expressed in terms of these operators as H^cs = J2j "^^j^j + ctes. The commutativity 
condition of the Ri operators follows from: i) the classical Yang-Baxter equation 



Kj' + [ri,j,ri^k\ + [ri,k,rj^k] = (12) 

where r^j = j^z:^ is the classical Yang-Baxter r-matrix, and ii) the equation \fi +'t^j,'fi^j\ = 
0. Unfortunately CRS, unaware of the Richardson's exact solution, gave not the eigenvalues 
Tj of the conserved quantities which were found in reference [10| using CFT methods, 



yU-1 J = l{7^«) 



8 Gaudin's model (1976) 



Inspired in part by Richardson's work on BCS, Gaudin proposed in 1976 a class of spin 



models based on a set of commuting Hamiltonians |11] 



^ t,:.t 



H^= J2 {i = l,...,n) (14) 



which can be diagonalized by an ansatz similar to (^) with the parameters satisfying 
the equations 



Nr. O , 



and eigenvalues 
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N 

E 



2ei 



(16) 



Compairing eqs. ( p^ and (^) it is obvious that Hi = —liuig^ooRi/gd, and similarly 
hi = —liuig^oori/gd. On the other hand, in limit g ^ oo, Richardson's eqs.(0) become 
Gaudin's eqs.(|T5|). The model defined by the Hamiltonians (^) is SU{2) invariant and it 
is known as the rational case. Gaudin also proposed a trigonometric model, which breaks 



SU{2) down to U{1) and an elliptic model which breaks this U{1) to Z2 |11]. Gaudin's 
trigonometric version have been generalized, to include a i^-term a la BCS, in references 



[12, 13 1 . There are also bosonic pairing Hamiltonians satisfying the same type of equations 



as in the fermionic case [14, 15[. Generalizations of the SU{2) BCS-like models to other 



Lie groups G |1£] and supergroups [17| have been worked out. 



9 Conformal Field Theory Picture (2000) 



In reference |1C] it was proposed a CFT interpretation and derivation of the exact's solu- 
tion of the BCS model. This was based on several observations: i) analytic structure of 
the Richardson wave funtion similar to the one that arises in the computation of confor- 
mal blocks, ii) common origin of the integrability of the BCS model and the Knizhnik- 
Zamolodchikov (KZ) equations [18[, namely the classical YB equations, and iii) similarity 
between the electrostatic analogue model of the Richardson's eqs. and the Coulomb gas 
representation of the Wess-Zumino-Witten (WZW) model [f^[ . 

The electrostatic picture was already noticed by Gaudin [20[ Richardson [^], who 
observed that the equations and ( p!5| ) are nothing but the equilibrium conditions for a 
set of N point-like charges, with charge Q = 2, located at the positions in the complex 
plane, subject to their mutual repulsion and the attraction of charges, with charge 
Q = —1, located at the positions Zi = 2ei, plus a constant electric field generated by a 
linear charge at infinity with density pL = —l/{7rgd). In the Gaudin's model the latter 
term is absent. The holomorphic piece, U, of the total 2D-electrostatic potential is given 

by 



U 



=1 l^u= 



:i ln(2;i 



Zj) 



u 



) + ^(-Ep=i^* + 2Elin.) 



(17) 



It is easy to verify that {dU / du^)^^^^^ = reproduces eqs.(|^ and (|T5|). Moreover the 
eigenvalues and hi are proportional to the forces {dU / dzi)^^^^^ exerted on the Q. fixed 
charges Zi. This analogue model was in fact used by Gaudin and Richardson to derive the 
standard BCS equations for the SC gap, chemical potential, total energy and occupation 
numbers of the energy levels, in the asymptotic limit when N ^ 00 and d ~ 1/A^ with 



N/n, kept fixed |i 
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From a CFT viewpoint one can recognized U as arising from the following chiral cor- 
relator |0| 



Q N 

= {W, n y-o^oi^^) n ^2aoK)) (18) 
i=l v=\ 

where K^(-z) and Wg are the following vertex operators 

= e*"'^(^), Wg = &-^^(^ ^dzLp{z)^ (19) 

constructed from a chiral boson </>(z) with background charge 2ao. In (^) we have ne- 
glected an operator at infinity needed to neutralize the overall background charge. We see 
from ( p^ ) that the Coulomb gas charges are equal to a = — ao for the Zj's, and a = 2ao for 
the iij/'s, in agreement with the Gaudin- Richardson electrostatic model, up to the overall 
ao factor. The line of charge at infinity is represented by the vertex operator Wg. 

The other ingredient of the CFT construction is the so called /? — 7 system formed by 
two boson fields fi{z) and 7(z) which have a correlator {(3{z)^{w)) = l/{z — w) |jl9|. Using 
this formula and the Wick theorem one can write the Richardson's wave function (P) as 
follows 

= (n.=i7^-"n^.) nli/?(e.)) (20) 

where mj = 1/2 ( resp. —1/2) if the corresponding energy level Zi is empty ( resp. occupied) 
by an electron pair. We have also neglected in ( |20| ) some extra /3— 7 fields placed at infinity. 
This suggest to associate every energy level to a primary field Vj^miz) of the SU{2) WZW 
model with total spin j and third component m, which have the Coulomb gas realization 
$^(2:) = ^^~"^{z) exp(—2ijao4>{z)). The total spin j is given by half the maximum number 
of electron pairs that can occupy a single energy level Zi, and so, in the actual case it is 
given by j = 1/2. Based on the eqs.(|l8|) and (pO|) it is natural to consider the following 
perturbed WZW conformal block (PWZW) 

i^^^^^{z) = {Wg^k{zi)---^L{^n){ duiS{ui)..J duNS{uN)) (21) 

where S{z) = P{z)exp{2iao(j){z)) are the screening operators needed to balance the charge 
and Cj/ are their contours of integration. Using the results explained above, one can write 
the PWZW conformal block as 

C^^^(z) = dm... duM e-"?. ^(^'") (z, u) (22) 

Hence, in the limit uq ^ oc these integrals are dominated by the saddle point configura- 
tions, i.e. {dU /duy)^ = 0, where the positions = satisfy the Richardson eqs., 



7 



and becomes proportional to V'm- Curiously enough, the proportionality factor C 

in eq.(^) follows from the gaussian integration around the saddle points, while the 1/\/nI 
factor arises naturally in the Coulomb gas approach. Given the relation k + 2 = 1/(2q!q), 
where k is the level of the Kac-Moody algebra SU{2)k, it follows that the limit ao — > oo 
corresponds to a singular limit in the representation theory of SU{2)k- Finally the PWZW 
conformal blocks (^) satisfy the perturbed KZ eqs. 

-^-Y Al^ I V;^^^^(z) = (23) 
2 dz, 2gd f^. z^-zjj^ ^ ' ^ ' 

In the limit when + 2 — > one can easily derive from ([2^ ) that ip^ is indeed an eigenstate 
of the operators Ri with eigenvalues ri = ^ gd dU/dzi (eq. (|l3|)). In this manner the 
Richardson's solution and the integrals of motion of BCS get unified in the framework of 
perturbed CFT. 

10 Gaudin's, BCS and Integrable Vertex models (1993- 
2001) 



The CFT interpretation of the BCS model, explained in the previous section, turns out 



to be closely related to the works of Babujian |22|, and Babujian, Flume [23| who in 
1993 rederived Gaudin's exact solution using the so called off-shell algebraic Bethe ansatz 
(OSBA). These authors also pointed out that Gaudin's eigenstates can be used to build the 
conformal blocks of the WZW models, along the same lines as was shown in the previous 



section. The work in reference |1C] was done with no knowledge of the papers |22, 23 1, 
which may then be regarded as another example of lack of information that these topics 
have unfortunately suffered in the past. 



The similarity between the works ||TQ] and [22, 53] suggested that the Richardson's 



solution of the reduced BCS model should also be derivable using the OSBA method. This 



was done by Amico, Falci and Fazio and later on clarified in references 26, 27 1, 
where the BCS coupling constant parametrizes a boundary operator that appears in the 
transfer matrix of the inhomogenous vertex model, whose semi-classical limit gives rise to 
the CRS conserved quantities. 

The OSBA approach starts from an inhomogenous vertex model whose transfer matrix 
is given by 



T(A;zi,...,Zf,) =tro(Koi?°''(A-zc)...i?°^(A-zi)) (24) 

where 

B^^ (A -Zj)=I^®h + r ® (25) 

ry - 2(A - Zk) 
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is the i?-matrix acting on the space (8> k ( are the spin matrices associated to the spin 
Sj irrep of SU (2) ) and 



A'„=(^'" 1] (26) 
is a boundary matrix acting on the auxiUary space 0. Defining the Bethe ansatz state 



\N) = B{Xi)...B{Xn)\ ],...,]) 



(27) 



one can follow two approaches: 

i) On shell approach: Impose that |A^) is an eigenstate of the transfer matrix T 
which is guaranteed by the Bethe ansatzs eqs. 



n 



An 



ry/2 + rjSi 



K- Zi- 7?/2 - r]Si 



N 

n 



Ag - A/3 + 7? 

Aq - A/3 - ?7 



(28) 



which in the limit r/ — > yield the Richardson eqs. 



Aq, 



Zi 



(29) 



The operators Ri (|11| ) appear at order r/ in the power expansion of the transfer matrix 

(H- 

ii) Off Shell approach: Let the state |A^) to be not an eigenstate of T but instead 
use it for finding solutions to the PKZ eq.(^). Then in the limit when k — > —2 one gets 
the Richardson's eqs. as saddle points conditions, as was shown in the previous section. 



For the Gaudin's model, this observation was made by Reshetikhin and Varchenko |28]. 
The derivation of Richardson's solution from an integrable vertex model has also al- 

] using the "determinantal" 



lowed the computation of correlators and form factors [29 



techniques developed in [30, 31 1. This works generalize the old results by Richardson and 
Gaudin's, concerning the norm of the eigenstates and the occupation numbers 



From a more mathematical viewpoint the CFT approach to the Gaudin model has 
a counterpart in the representation theory of Kac-Moody algebras at the critical level, 



as was shown by Feigin, Prenkel and Reshetikhin [32|. These authors remark that their 
construction fits into the program of geometric Langlands correspondence proposed by 
Beilinson and Drinfeld. There are also connections with the Gaudin-Calogero model and 
the Hitchin's systems p3| , An interesting mathematical problem is to investigate the 
meaning of the Richardson eqs. in connection with the representation theory of Kac-Moody 
algebras at the critical level. 
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11 BCS and Chern-Simons theory (2001) 



In this section we shall briefly mention the last step in the understanding of the integrability 
properties of the BCS model. In the Gaudin's case we can draw the following chain 
of relations: Gaudin's magnets Integrable vertex models — > WZW — > Chern-Simons 
theory. The last step refers to the well known connection between the 3D-CS theory and 
the 2D-WZW model [35|. Hence we may ask what is the origin of the BCS model in the 
CS theory? 

In reference [16| it is shown that the field theoretical origin of BCS can be traced back 
to a SU{2) CS theory interaction with a one-dimensional distribution of colored matter 
which breaks both gauge and conformal invariance. This connection is quite remarkable 
because Chern-Simons theory has been advocated to be mainly connected with effective 
descriptions of fractional quantum Hall effect and high superconductivity, but never 
with standard superconductivity. The Chern-Simons theory is not defined in the physical 
space, which might be two or three-dimensional, but rather in the complex energy plane 
which is always two-dimensional. This explains why this field theoretical connection of BCS 
theory remained unveiled for so long time. The connection of Chern-Simons theory with 
BCS model can be understood in a more general framework when one considers a scaling 
limit of the twisted Chern-Simons theory defined on a torus, that is, the twisted elliptic 
Chern-Simons theory. On a torus the KZ equations [18| are replaced by the Knizhnik- 
Zamolodchikov-Bernard (KZB) equations [^], which depend on the coordinates Zi of the 
punctures, the moduli of the torus r and a set of parameters Uj characterizing the toroidal 
flat gauge connections. The later parameters uj define the twisted boundary conditions 
for the WZW fields on the torus. 

The main result of [^] is to show that, for a generic simple simply connected, compact 
Lie group G, the Richardson equations, the CRS conserved quantities and their eigenvalues 
arise from the KZB connection and their associated horizontal sections. This is done in a 
limit where the torus degenerates into the cylinder and then into the complex plane. In this 
limiting procedure the generalized BCS coupling constants appear as conjugate variables 
of the parameters uj, when this parameters go to infinity. This gives the G-based BCS 
models a suggestive geometrical and group theoretical meaning. 

In analogy with Gaudin's model we can complete the chain of relations in the BCS case: 
Richardson solution Integrable vertex models with boundary operators — > perturbed 
WZW — > twisted elliptic Chern-Simons. These connections suggest further generalizations 
both of the CS theory and the CFT, which deserve further investigation. 



12 Summary 



In this article we have summarized some of the works concerning the integrability and 
conformal properties of the exactly solvable reduced BCS model and the closely related 
Gaudin's model. 
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